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Introduction 


All  tli©  properties  wlitch  describe  the  state  of  an  explosiv 
substance  Immediately  after  the  passage  of  a  steady-state  one- 
dimensional  detonation  wave  oan  be  oaloulated  with  considerable 
accuracy  by  the  hydrodynamio-thermodynamio  theory.  This  theory 
determines  the  five  unknown  properties 


fl)  Detonation  velocity  D 

(2)  Gas  velocity  U 

(3)  Detonation  temperature  Tg 

(4)  Detonation  pressure  3?g 

(5)  Detonation  spec if to  volume  vg 


in  terms  of  the  known  thermal  properties  of  the  initial  state,  with 
the  aid  of  the  five  equations 


) 

2) 

3) 

) 

5) 


Conservation  of  mass 
Conservation  of  momentum 
Conservation  of  energy 
liquation  of  state  of  the  substance 
The  Chapman -Jouguet  condition, 

V2Li  s  (^s) 

Vn-V,  ~  \dVr,  ii 

iJ  JL  *  *J 


These  steady-state  properties  depond  only  upon  the  equilibrium 
conditions  in  front  of  the  wave  and  behind  the  wave,  and  are 
completely  independent  of  the  detailed  structure  of  the  detonation 
wave  and  of  properties  which  affect  only  that  structure  — 
properties  like  the  viscosity,  tho  thermal  conductivity,  or  the 
specific  rate  of  the  chemical  reaction. 

Various  researches  on  the  necessary  high-ter.iperature  thermal 
data  and  on  the  appropriate  equation  of  state  for  highly  compressed 
gases  (of.  the  summaries  by  Kistiakowsky  and  Wilson,  and  by  Cook) 
and  on  the  validity  of  the  Chapman-Jouguet  condition  fvon  Neumann) 
permit  calculations  of  tho  detonation  properties  of  any  chosen 
explosive  to  be  made  with  considerable  assurance.  In  the  present 
investigation  we  shall  assume  the  correctness  of  the  hydrodynamio- 
thermodynamio  theory  for  steady-state  one -dimensional  waves,  and 
shall  use  its  results  whenever  necessary. 


The  agreement  between  theory  and  experiment  becomes  poor 
for  some  explosives  under  certain  conditions  —  in  particular,  for 
explosives  of  low  energy,  loosely  confined  or  in  narrow  charges, 
in  the  cast  state  or  large  grains,  or  insufficiently  boosted. 
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Under  these  conditions  the  explosive  nay  fa)  refuse  to 
or,  fall  to  attain  the  full  theoretioal  detonation  velocity, 
problem  is  to  decide  whether  a  particular  initial  impulse  will 
attenuate  and  eventually  die  away,  or  whether  it  will  grow  into 
a  detonation.  In  other  words,  a  solution  of  this  problem  would 
cover  the  entire  subject  of  "sensitivity"  of  explosives  (i.e., 
sensitivity  to  thermal  initiation.  The  problem  of  sensitivity  to 
shook  involves  the  additional  problem  of  the  conversion  of 
meohanloal  energy  into  heat). 

Or  the  explosive  may  (b)  after  detonation  has  been  initiated,  fall 
to  maintain  it.  After  the  explosion,  part  of  the  explosive  remains 
unchanged.  The  solution  of  this  problem  requires  a  study  of  the 
rates  of  the  dissipative  processes  whloh  cause  the  ohemioal  reaotlon 
(and  the  detonation  wave)  to  oome  to  a  stop. 

The  present  Investigation  has  been  undertaken  to  study  the 
theory  of  the  causes  of  these  two  kinds  of  failure.  In  the  oourse 
of  It  we  shall  study  the  distribution  of  temperature  and  pressure 
within  a  detonation  wave,  and  the  effect  of  these  variables  upon 
the  rate  of  the  ohemioal  decomposition  reaotlon. 


The  report  will  be  divided  into  four  parts; 

Part  I.  Ohemioal  Reactions  in  an  iflcploslve  Substance. 

The  effect  of  temperature,  pressure,  and  the  nature  of  the 
reacting  phase  will  be  discussed. 

Part  II.  Detailed  Structure  of  the  Steady-State  Detonation  wave. 

Wumerioal  integration  will  be  employed  to  obtain  the 
distribution  of  temperature  and  pressure  in  shook  waves  with  and 
without  ohemioal  reaction. 

Part  III.  Detonation  Wave  with  Radial  Losses. 


Part  IV.  Detonation  Wave  in  the  Non-Steady  State. 


X00&*’ 
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'Jlne-Teuperature  Relation  In  an  Explosive  Reaction. 
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It  will  be  useful  to  have  a  preliminary  idea  of  the  times 
required  to  attain  various  temperatures  in  the  explosive  reaction. 

'*  ?o  obtain  this  preliminary  estimate,  we  ean  consider  a  volume  of 
reactant  all  of  whioh  is  brought  to  the  temperature  T0  at  time  t=0. 

We  then  solve  the  conservation- of -energy  equation  with  the  neglect 
of  all  transport  terms,  obtaining  the  temperature  T  as  a  function 
of  time  £.  (it  is,  of  course,  realised  that  the  neglected  transport 
terms  may  become  important  in  an  actual  wave. ) 

Consider  the  enrgy  change  within  this  reaction  mixture  from 
time  t  •  0  to  time  t  s  £.  Then 


f3  -  aQ)+  Q  -  CL (T  -  t  )  =  0 


where  T  is  the  temperature  at  time  jb,  ...... 

is  the  mean  heat  capacity, 
is  the  heat  evolved  by  the  reaction  to  time  t, 

(a-aa)  inoludes  all  other  energy  terms:  radiation,,  conduction', 
kinotio  energy,  adiabatic  compression.  It  will  in 
general  be  a  function  of  both  time  and  distance.  - 

To  the  approximation  we  now  require,  (S-Ep)  may  be  set  equal  to' 
zero,  giving 


T  -  In  5  .Q/S_ 


*»  ■ 


By  neglecting  all  transport  terms  we  have  eliminated. distance-  from 
the  equation.  In  other  words,  it  is  assumed  that  all  the  heat  of 
the  reaction  goes  to  raiBe  the  temperature  of  the  reaction  mixture* 
At  infinite  time  t  =  oo  f  the  temperature  will  have  attained 
a  terminal  value 


Too"  T0  = 


. (3) 


where  is  the  heat  evolved  in  the  complete  reaction. 
•  The  temperature  variation  is  now  given  by 


i 


-  4  - 


To  proceed  farther,  We  must  evaluate  dQ/dt  in  terms  of 
kinetios  of  the  reaction.  ;*or  a  first  order  reaction 


'%4i 


■  i.H 

at 

n  .  -fze^dt 
n  "  0  c 


(6) 

(6) 


where  n  is  the  amount  of  unchanged  reactant  remaining, 

Z  is  the  probability  factor  for  the  reaction, 

A  is  aH/R,  the  heat  of  aotivation  for  the  reaotion  divided 
by  the  gas  oonstant. 

The  heat  Q  is  then 


and  its  time  derivative  is 


d£  . 

at 


^oo  « 


(0) 


On  combining  £q.(7)  and  ©1.(0), 


•  teoo-  <!)  Z  e"AA  . (9) 

The  temperature  variation  is  therefore  given  by 

ff  =  (Too-  S)  Z  e-V®  . (10) 


which  is  the  basic  differential  equation  to  be  solved.  It  is 
integrated  without  difficulty  to  give 


where  the  symbol  Ei  represents  the  exponential  integral 
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for  whloh  tabulated  values  are  available.  For  all  valued 
temperature  except  those  in  the  immediate  neighborhood  of 

the  exponential  integral  may  be  replaced  by  its  approximate  value 


El  x 


s  giving  as  the  final  result 


t 


1_ 

AZ 


2oo-T 


0 


_rf_  a/i 

Too-T  °  . 


(12) 


From  Eq.  (18 )  —  orf  in  the  neighborhood  of  Tqj  ,  from 

Eq.  (11)  —  the  time- temperature  curve  fer  the  entire  course  of 
the  explosive  reaction  oan  bo  calculated  readily.  As  an  illus¬ 
trative  example,  this  has  been  done  using  the  following  numerical 
valuos* 


T0  B  1000°2 
i’oo  =  6000°K 
Z  =  1016  boo-1 
k  -  60,000/1.986 

and  the  results  are  tabulated  in  fable  I.  The  same  results  are 
graphed  in  Fig.  1,  and  to  a  greatly  enlarged  scale  in  Fig.  £. 


*The  value  of  Z  was  taken  to  be  kg/h  where  k  is  Boltzmann's 
constant,  T  is  3000°IC,  and  h  is  Planck's  constant,  multiplied  into 

e^/R  with  =  5  oal  dog"1.  These  are  reasonable  values,  for  the 
entropy  ohange  of  activation  in  the  highly  condensed  system  is 
probably  small.  ± 

The  value  of  aHr  was  arbitrarily  ohosen  as  60  kcal.,  or 
about  the  strength  oFa  carbon-carbon  bond.  Experimental  values 
from  low-temperature  measurements  have  been  reported  both  higher 
and  lower  than  this  value. 
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hiapussion 

fa)  The  calculates  curves  will  be  observer  to  reproduce 
the  qualitative  matures  Ox’  an  explosive  reaction*  a  long  induction 
period  is  Tolloweu  by  an  extremely  rapid  rise  o2  temperature,  which 
does  not  cease  until  the  reaction  is  substantially  complete. 


(b )  Y/hen  the  reaction  lias  proceeded  to  any  appreciable 
extent,  the  second  term  on  the  ripht  or  Oq.  (12)  may  be  neglected 
in  comparison  with  the  rirst,  ;iviny; 


-  1 


•00"  i?0 


(!£>) 


The  induction  period  therefore  depends  exponentially  on  the 
reciprocal  oT  the  initial  absolute  temperature.*  This  is 
precisoly  the  rosult  round  experimentally  by  vaiioiiB  experimenters 
( bn  dr  gov,  Garner,  Katies,  Harvey,  Iiuraour).  The  relation  is 
usually  written  in  bhe  Torm 


log  t  =  4  +  B  . (14) 

■4. 


*In  the  illustrative  example,  the  induction  period  ror 
.*?0  =  300°i:.  is  10*  *  years.  Tho  substance  would  be  stable  at 
room  temperature 2 


TONFlDENf'W- 
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Jills  agreement  with  experiment  leaus  us  tc  expect  tha1 
foregoing  time -temperature  analysis  ie  not  far  from  oof 


'.$!/$ fi. 


wr/A 


(o )  The  numerical  values  oi‘  Jigs.  1  and  2  are  probably 
sufficiently  aoourate  to  provide  an  order- of -magnitude  estimate. 

It  will  be  observed  tliat  the  last  75<b  of  the  reaction  requires 
only  about  10"11  seconds.  During  this  time  a  molecule  In  the 
compressed  state  will  undergo  some  10  to  100  collisions,  whloli  will 
serve  to  equalize  any  local  temperature  fluctuations. 

hVen  with  an  Initial  temperature  of  1000°K.,  (which  ie  a 
temper  -cure  tliat  may  reasonably  be  oirpectoa  to  occur  looally  from 
frictional  heating)  the  Induction  period  is  only  of  the  order  of 
10~5  seconds.  It  is  a  fact  that  time-lags  Oi  ignition  oi  sensitive 
oxpl os ives,  du  determined  from  oscillograph  records,  aio  about 


1G"5  seconds  fcuoted  by  Cook,  194ti). 

If  olio  reaction  does  prooeed  at  velocities  of  ohe  order  of 
those  estimated  here,  there  remains  no  objection  to  tho  "thermal" 
^ypo^hosis  or  reaction;  the  only  previous  objection  liad  been  tliat 
the  reaction  wo^ld  prooeed  too  slowly  (e.g.,  one  estimate  gavo  a 
time-lag  of  10" J  seconds). 

\'»*e  may  therefore  make  subsequent  calculations  on  tho 
assumption  that  each  individual  molecule  in  the  explosion  undergoec 
an  ordinary  thermal  reaction,  and  that  its  rate  of  rotation  can  bo 
treated  by  the  usual  methods  of  ohomioal  kinetics.  In  particular, 
it  will  not  be  necessary  to  assume  that  the  reaction  is  initiated 
by  any  special  mechanism  Involving  impact  or  shock  waves,  as  has 
sometimes  been  assumed. 


*Semenov  lias  obtained  an  equation  also  of  the  form  of  :sq.(14) 
by  assuming  .that  the  reaction  produces  fragments  which  can  later 
initiate  chains  —  the  so-called  "degenerate  branching. 11  It  is 
undeistondable  tliat  the  equations  should  be  formally  similar 
because  in  each  process  something  produced  by  the  reaction  (heat 
or  chain-starters,  respectively)  serves  to  accelerate  the  reaction. 
Che  "degenerate  branching"  hypothesis  seems  less  likely  than  the 
present  "thermal"  hypothesis  for  detonation  in  solid  explosives, 
since  tne  detonation  process  has  not  been  shown,  to  be  sensitized 
by  added  cnain-initiators . 
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Ex  feet  of  Pressure  on  an.g.lfirplosive  Reaotlon 


We  now  proceed  to  investigate  the  consequences  of  our  r// 
assumption  that  the  elementary  rate  process  in  an  explosive  react r 
may  be  treated  by  the  usual  methods  of  roaotlon  kinetics.  The 
t  absolute  rate  of  a  chemical  reaction  la  given  by 

^  =  S  . (IB) 


where  kr  Is  the  specific  reaction  rate  constant, 
k  is  Boltzmann's  constant, 

T  is  the  absolute  temper  at  me, 
h  is  Planck's  constant, 

is  the  standard  free  energy  change  for  the  formation  or 
a  mole  of  activated  complex  from  the  reactants, 

R  is  the  gas  constant  in  the  same  units  as  ?/j . 


To  investigate  tho  effect  of  pressure  P  on  the  specific 
reaction  rate,  v/e  wish  to  calculate  the  increase  in  free  energy 
when  the  system  goes  from  its  normal  to  its  activated  state  at 
the  pressure  P.  The  thermodynamic  relation  is 


i  DiA  - 
VSFb  " 


which  lute'. ’rated  bocoueu 


*  1  (  V 

xi’p-j,  =  a.'? s0  *■  Uv*ar 


=  i? 


This  is  the  increase  in  free  energy  of  activation  which  results 
from  the  work  done  in  the  volume  change  of  activation  Ay  against 
the  external  pressure  P.  It  seems  safe  go  assume  that  in  the 
decomposition  reactions’  found  in  explosions  the  activated  state 
will  be  more  expanded  in  volume  than  the  reactants.  Thus  AV*  is 
always  positive,  an.:  a  reasonable  estimate  of  its  magnitude'  is 
10  njiW  mole1"1  (of.  b tears  and  Eyring).  If  we  take  the  pressure 
in  an  explosion  to  be  50,000  atmospheres,  the  increase  in  free 
energy  is 

JL.  r 

p =  50o, OUO  cmJ  atmospheres 


=  IE  kilocalories 
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Since  the  volume  change  Oi  activation  of  a  dilute  gas  is  z 
this  freo  energy  is  also  the  di  if or once  between  the  free  ener 
oi  activation  in  the  gas  phase  an'.  In  the  liquid  (or  highly 
compressed  gas)  ph^so; 


liquid  -  =  P^S  12  k°al  . (19) 

At  a  temperature  of  5000°K.,  this  dlfioronco  in  froe  energy  of 
activation  Is  enough  to  slow  the  reaction  rate  down  by  a  faotor 

of  e1*2  s  5.  fhe  effeot  of  pressure  on  the  specific  rate  at  high 
temperatures  is  therefore  so  small  that  it  may  safely  be  ignored 
in  all  but  the  most  highly  refined  oaloulat i ons • 

Sinoe  the  sign  of  the  volume  ohange  of  activation  will 
most  probably  be  positive,  the  effect  of "high  pressure  will  always 
be  to  slow  the  reaction  down  if  the  reaction  is  unimoleoular . 

An  increase  of  pressure  would  novor  be  expected  to  have  the  effect 
of  markedly  i nor easing  the  specific  tfeaction  rate,  as  has  sometimes 
been  pr  op  os  e3T  * 


Homogeneous  Reaotion  in  the  Solid  State 


Suppose  the  activated  state  ie  to  be  formed  in  the  interior 
.  o*  a  solid*  Since  the  reaction  has  been  supposed  to  reqtire  a 
*  volume  increase  on  activation,  there  will  be  a  corresponding 

increase  In  the  free  energy  of  activation,  becauso  of  the  additional 
work  of  compression  on  the  surrounding  solid.  One  way  of  estimating 
the  increase  in  free  energy  is  to  suppose  that  the  surrounding 
solid  contracts  homogeneously.  Again  making  use  of  the  thermo¬ 
dynamic  relation 


(16) 


we  integrate  to  obtain 


*  whore  ^  is  tho  compressibility  of  the  solid,  fhis  is  the  increase 
in  free  onergy  of  activation  which  results  from  the  work  done  in 
the  deformation  of  the  solid  lattice,  fhe  compressibility  in 

.»  kg"1  cm2  of  fats  and  related  compounds  is  about  50  10"6,  and  of 
granite  is  about  £  10~6.  fhus  the  increase  in  froe  energy  is 


\ 


/ONFIDENI'. 


Cl 


\ 


■  10 


106  kg  om 
6  5 

10  om  atmospheres 


24  kilocalories 


"ince  a  liquid  ie  sufficiently  disordered  that  its  lattice  will 
iot  ho  compressed  like  that  of  a  solid,  this  free  energy  is  also 
*  the  difference  between  the  free  energies  of  activation  in  the 
.-solid  phase  and  in  the  liquid  (or  compressed  gas)  phase: 


^6 olid  -  ^liquid 


>  24  teal 


(21 ) 


While  at  room  temperature  this  difference  in  free  energy 
would  cause  the  rate  to  he  slov/er  in  the  solid  state  by  a  factor 

of  e40«  10«  and  so  cause  the  solid  to  be  enormously  more  stable 
than  the  liquid  or  gas9  at  a  temperature  of  4000°£.  the  reaotion 
in  the  solid  will  only  be  slower  by  a  faotor  of  e^  or  about  20. 
fhus  if  the  temperature  is  high  enough,  the  solid  will  decompose 
at  a  rate  comparable  to  that  of  the  liquid  or  vapor. 

Discuss  ion 

flie  parallel  between  the  decomposition  reactions  of 
explosives  in  the  gas,  liquid,  and  solid  states  and  diffusion 
processes  in  those  same  states  is  complete. 

In  the  dilute  gas  phase,  the  diffusion  process  goes  without 
activation  energy  and  does  not  depend  on  pressure.  I'he  volume 
change  necessary  in  the  decomposition  reaction  can  take  place  at 
no  expense  of  free  energy,  since  each  molecule  has  plenty  of  space 
already. 

In  the  liquid  phase,  a  diffusing  molecule  requires  a  certain 
amount  of  additional  volume—  or,  as  it  lias  been  called,  a  "hole? 
J?he  formal} ion  of  this  hole  requires  a  certain  amount  of  free 
energy  (a  fraction  of  the  heat  of  vaporization)  and,  in  addition, 
the  free  energy  required  to  carry  out  the  expansion  against  any 
external  pressure  P.  In  the  decomposition  reaction,  free  energy 
must  similarly  be  supplied  for  the  volume  change  against  the 
external  pressure  P, 

In  the  solicT  phase,  a  diffusing  molecule  also  requires  extra 
volume  to  permit  its  flow.  But  where  in  a  liquid  such  a  hole 
could  be  introduced  without  appreciably  disturbing  the  order  of  the 
entire  liquid  (since  it  was  already  disordered),  in  a  solid  the 
introduction  of  such  a  hole  produces  a  displacement  of  a 
considerable  number  of  molecules  in  the  lattice  from  their  previous 
positions.  This  requires  more  energy,  consequently  diffusion  in 
the  solid  state  is  much  more  diffioult  than  that  in  the  liquid 
state.  j?he  same  is  true  of  decomposition  in  the  solid  state.  It 
is  tho  requirement  that  the  lattice  retains  its  order  —  that  a 
great  many  molecules  in  the  lattioe  must  move  cooperatively  when 
the  lattioe  expands  —  that  suppresses  both  diffusion  and  ohemioal 
reaction. 
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Tho  Transition  Between  Reap t ion  Meohanisms 


Both  from  tho  considerations  of  the  previous  eeotlou  which 
Indicate  that  reaotlon  should  bo  possible  in  the  solid  state,  and 
from  the  experimental  fact  that  it  is  possible  to  detonate  oast 
8 ticks  of  explosive,  we  are  lod  to  write  the  reaction 


u 


RSVBiUIU  >  PRCDUCjB 


(22) 


It  is  also  known  that  under  conditions  where  dissipation 
is  possible,  detonation  velocity  is  very  sensitive  to  partiolo 
size,  smaller  partiolos  (i.e.  more  total  surfaoe)  giving  velocities 
nearer  to  tho  theoretical,  The  conclusion  immediately  suggests 
itself  that  the  reaction  is  a  surfaoe  (or  "topoohemical" )  2/3  order 
reaction,  evaporation  probably  being  the  rate-determining  step. 

The  reaction  may  be  formulated 


\ 


RJilCTAUT 


kovap  v 
*  ^<5v ap 


VAPOR 


VAPOR  kg^e  > 


phosuoi’s 


(23) 


Wo  now  enquire  when  reaction  (22)  will  outstrip  the  system  of 
reactions  (23).  The  kinetics  of  reactions  (23)  may  be  obtained 
by  assuming  a  steady-state  concentration  of  vapor.  Then 


d  (VAPOR) 

dt  "  kevapteURAACifi)  -  k4vap(SUR^Clfl)  (V*j?OR) 

-kgae  (VaPOR)  =0  . (24) 


keVap  ( ^UR..*!<  «Cii) 


^evap  ( &  OB  ?AC  U ) 


-dfPRODOOMfi)  .  . 

— dt — rrrr..(2D) 


+  i 


O 


gas 


The  two  limiting  casos  are  of  racst  interest;  If  evaporation  is 
-  much  slowor  than  gas  reaction, 


d  (PRODUCTS) 
dt 


kovap(smi,W:i 


(26) 
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If  on  the  other  hand  the  sac  reaction  is  slowor, 

d (PRODUCE  )  kovap 

dt  *  'I^8  ^ovap 

Above  tho  temperature  at  which 


kgas  =  kevap (SURjIAG  Ifi ) 

the  rate-dot ormini ns  process  is  evaporation;  below  that  temperature 
it  is  the  gas  reaction.  In  order  to  find  tho  absolute  rates  of 
those  two  processes,  wo  estimate  their  heats  and  entropies  of 
activation  and  apply  the  absolute  rate  theory; 

4 

,  „  kf  -^electronic 

Ve  "  ~h  0  - I5T .  . (28) 


The  only  activation  onorgy  for  tho  C'5-6  reaction  is  that  nocoesary 
for  the  electronic  rearrangement  which  accompanies  tho  reaction, 

kova?fBURVAOJ)  =  lO-6  ~  e10/K  . (29) 


Tho  fraction  of  moleculoe  which  are  in  tho  surface  is  really 
variable,  since  it  depends  on  tho  2/3  power  of  the  amount  unreacted. 
The  fraction  hero  taken,  10“ 6 ,  is  an  estimate  which  might  vary 
from  case  to  caso.  Condensation  is  assumed  to  go  without  activation 
onorgy,  but  requires  tho  loss  of  the  entropy  of  vaporization,  which 
is  estimated  to  be  10  cal  deg"-*-  for  the  oompr eased  gas. 

VThen  the  two  rates  are  equal, 

* 

=  ^electronic _ 

K  In  6  +  &S°yap  . tB0)  ■ 


m  30. 000 
*  2V.4  4  10 


=  G00°II. 


The  electronic  activation  energy  is  estimated  to  be  30  kcal.  The 
final  result  is,  of  course,  subject  to  considerable  error. 


The  rate  of  the  homogeneous  reaction  is 

v  .  v--,  -  f  ^electronic d^o^ansion) 

oOlid  -  e - 55 - 


(Si) 


confidential 


where  the  fraction  of  surface  molooules  is  taken  as  before,  and 
evaporation  is  assumed  to  require  tho  full  heat  or  evaporation 
and  to  proceed  without  entropy  of  activation. 

Reactions  (22 )  and  (26)  nov/  have  tho  Eamo  rate  at  a 
temperature 


f  2  ~  e3-ec  +  6  jloacpans  i on 

.  R  In  6 


d^eyap _ 


(ob) 


530,  ooo  4  30,000  -_io,opo 

27.4 


=  iaoo°£. 


'■’he  ii*ee  energy  necessary  ror  expansion  is  estimated  to  be  bO  kcal, 
llie  heat  of  vaporization  10  kcal.  This  final  result  i.lso  is 
subject  to  considerable  error. 


%  V/e  may  now  assemble  the  rosults  o^.  the  abo\e  calculations. 

lat )Solid  reaction 

•  Above  ...  temperature  or  about  2000°ic. ,  reaction  v/ill  by 

preference  occui  in  the  solid  phase,  homogeneously •  The  late  will 
obey  Ifiq.  (531).  At  hi  :li  tempera lures #  the  hi  ;i  ac  uivat  ion  energy 
is  not  a  particularly  serious  handicap  to  the  reaction  (see  Pigs. 

53  and  4),  In  fact,  at  high  temperatures  the  reaction  is  not 
markedly  temperature -dependent  (  'i  ;.  4).  however  at  lov/er 
temperatures  the  reaction  becomes  highly  temperature-dependent  and 
drops  out  completely.  Vais  means  that  if  the  reaction  iB  to  go 
in  the  solid,  tho  temperature  must  be  kept  high,  The  mean  length 
of  the  reaction  zone  is  10&§  upward  (see  pigs.  53  i*ncl  4). 

(b )  hVap oration 

Prom  about  2000°!:.  down  to  about  G00°i:.,  the  preferred  rate- 
determining  step  is  evaporation.  It  obeys  ,jq.  (532)  and  is  only 
mildly  tenporatur e-dependent,  iji^ce  it  is  proportional  to  the 
amount  of  surface,  this  reaction  rate  will  be  sensitive  to  particle 
c  size,  fho  reaction  zone  is  of  the  order  of  1  to  100  centimeters, 
so  this  reaction  rate  v/ill  also  be  sensitive  to  the  size  and  form 
of  the  explosive  charge.  If  the  roaction  is  to  go  rapidly  by  this 
nochan  ism,  the  particle  size  must  be  kept  small  or  the  charge  size 
*  lt-rge  • 

->Q  )  Pas  reaction 

This  is  the  only  reaction  below  about  G00°i:.  It  is  so 
temperature-dependent  c*nd  the  ieaotion  zone  so  long  that  there  is 
little  chance  of  it  going  at  all.  In  fact  Belyaev 
explosion  initiation  does  not  occur  until  the  boil 


finds  that 
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Three  types  oi’  reaction  in  an  explosion.  Rate  plotted 
logarithmically  against  temperature. 
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logarithmically  against  reciprocal  absolute  temperature. 
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PART  II.  JHTAILiSD  STRUCTURE  01?  THE  6 T^ADY-b JSATli  1)  WAVE 


The  Role  of  Thermal  Conduntion 


It  will  be  of  int erect  to  eee  whether  thermal  conduction 
„  is  of  any  Importance  in  the  propagation  of  the  detonation  wave. 

To  obtain  an  approximate  answer  to  this  question,  we  take  a  semi- 
infinite  medium  of  thermal  conductivity-  k>  ,  impose  on  it  the 
boundary  condition  All  =  AT-.  at  x=0  for  all  time  £3*0,  ancl  follow 
the  temperature  witETn  tEe  medium  as  a  function  of  time  and  distanjo 
This  is  a  classical  problem  and  the  solution  is  well  known 
(of.  Carslaw,  Conduction  of  Heat  p.  46),  If  the  differential 
equation  for  t he  dl ffus ion  of  temperature  is 


(34) 


the  required  solution  is 


2  Erf 


where  AT  la  the  temperature  riso  at  any  x,  £p 

aTi  is  the  temperature  rise  at  the  boundary  x=0, 

K,  is  the  thermal  conductivity, 

V  is  the  volume  (per  gram,  mole,  or  raoleoule), 

Cv  is  the  heat  capacity  (for  the  same  quantity), 

~~  .  'X  2  . 

-  -  rr—  \  -y  /2 

Erf  x  is  the  error  function  ~  j  e  dy 

v  4  'K  V 

for  vdiioh  values  ore  tabulated. 


(35) 


We  now  require  a  numerical  value  for  the  parameter  leV/Oy • 

.  straightforward  application  of  kinetic  theory  (of.  Jeans, 
Dynamical  Theory  of  Gases  )  gives 

*  fce  -  1/3  W/v  CyO  1-  (persistence  of  velocity  factor)  ••••(36) 

,  *  whore  n/v  is  the  number  of  molecules  per  unit  volume, 

7  is  the  moan  velocity, 

I.  Is  the  mean  free  path. 
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?or  a  liquid  this  becomes  at  onoe  (Kinoaid  and  gyring) 
^  _  (persistence  faotor)  .  ,.3  _ 

** - 3 - L  1/X  UA 


(37) 


where  u  is  the  velooity  of  sound, 

X  ie  the  into  rat  omio  distance. 


We  shall  use  this  expression  for  the  highly  compressed  gas  also. 

The  quantity  K>v/cv,  which  may  be  called  the  "thermal  diffusivity, n 


icy  -  (persistence  faotor)  \ 

G  3 

v 


(36) 


:*or  a  gas  with  spherical  molecules,  the  persistenoe-of-velooity 
factor  is  £.622  (Chapman  and  Cowling,  p.  £35).  With  the  high 
forward  velocities  existing  in  the  detonation  wave,  this  factor 
might  be  somewhat  larger.  In  any  event,  it  will  approximately 
cancel  the  factor  1/3  so  that 


(39) 


The  diffusion  equation  and  its  solution  for  the  problem  of 
the  semi-infinite  medium  now  become,  respectively 


&  =  UX 


(40) 


2  Erf 


\fz  UA 


With  the  numerical  values  u 


6  106  cm  sec 
3  10~8  cm 


(41) 


the  temperature  rise  has  been  computed  from  i3q.  (41)  for  various 
times  and  distances.  The  results  for  1  A,  10  A,  and  100  A  are 
plotted  in  ?ig.  5,  and  the  numerical  results  for  10  A  are  also 
given  in  Table  II. 
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The  range  of  effective  thermal  conduction  in  a  detonation. 
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Table  II 


Time,  seconds 


10' 

10' 

10' 

10' 

10' 

10' 


13 

12 

11 

10 

9 

0 


0.096 

•  698 

•  868 

•  958 

•  9867 

•  9958 
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If  the  detonation  wave  is  traveling  with  a  velocity 

D  =  8  10®  cm  sec"1,  the  temperature  at  100  S.  distance  would  not 
nave  been  raised  peroeptibly  by  thermal  conduction  in  the  time 
required  for  the  wave  to  travel  100  a*  Therefore  at  this  distanof. 
thermal  conduction  is  of  no  importance* 

At  10  A •  thermal  conduction  will  have  raised  the  temperature 
about  15$  of  its  total  value,  in  the  time  the  wave  would  take  to 
travel  this  far.  Therefore  at  this  distance  thermal  conduction 
is  of  a  certain  importance. 

At  1  A,  thermal  conduction  will  have  raised  the  temperature 
65 fo  of  its  total  value  and  is  therefore  of  great  importance. 


The  effect  of  conduction  estimated  in  this  way  is  probably 
an  upper  limit,  einoe  in  a  reacting  system  the  source  of 
Temperature  would  not  maintain  its  high  value  continuously,  nor 
would  the  velocity  of  sound  always  be  so  high  as  that  assumed  here. 
This  leads  us  to  the  conclusion  that  the  conduction  of  energy 
fsuch  as  that  given  out  by  a  reacting  grain  of  explosive)  in  the 
compressed  gas  is  comparatively  short-range. •  .at  most,  over  a  dis¬ 
tance  of  a  few  molecular  diameters.  Thus  the  temperature  within 
a  reacting  system  will  not  be  the  some  throughout.  It  seems  more 
likely  that  the  temperature  in  the  given  region  will  not  bo  very 
different  from  that  resulting  from  the  ohomioal  reaction  whioh  hiib 
taken  place  in  that  given  region.  Those  regions  in  whioh  react  lor. 
has  prooeodod  farther  v/ill  then  bo  at  a  higher  temperature. 

If  this  qualitative  result  is  correct,  the  temperature 
distribution  in  a  detonation  wave  with  ohemloal  reaction  would 
more  nearly  resemble  figure  6a  than  Figure  6b  which  is  usually 
assumed; 
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Thu  Dlfforont ial  aquations  for  a  Ono -Dimensional  steotly- 
State  Detonation  Wave 


The  propagation  ox  a  one-di mend o.nal  detonation  v/ave  is 
donor  theft  "by  olio  partial  differential  equations 


fc+»fc  =  -<$ 


for  tho  co  u  a  or  vat  ion  of  mass, 


+  T!&y  -  “  \p 


f oi  tho  conservation  of  momentum,  and 


^  ,  ujyj  -  i  (TJ  u2ur^'  4  i  ^  (l/O'j? 

+  _  "e  V  f  jfV% 


“Su  \  .  l  3 


. (43) 


for  tho  conn  oi  vat  ion  of  onorcy.  Ilorc  m  Lb  4/3  times  the  vie  cos  it: 
and  other  symbols  '^avo  already  ho  on  defined. 

If  tho  detonation  wave  is  moving  a  Ion:;  at  a  constant 
velocity  and  with  a  constant  shape  and  magnitude,  the  trans  format 1 
to  a  coordinate  oyster*  moving  with  tho  velocity  of  the  wave  purmi 
tho  conservation  equations  to  ho  integrated  at  once  (of.  Loolcor  o- 
ICiotiakowslcy  and  7ilson)  to  obtain  who  total  di  f  for  out  ial  oquatio 


>u  -  iiL 


I-lu^P  -  J  "  ya  dn/dy 
/ 

B  +uJ/ii  -  fa*4  -  x-'  -  ~  iu'/dy 


where  u  is  the  velocity  in  the  movin';  coordinate  system, 

2  is  tho  distance  in  the  moving  coordinate  uyavein, 

ii,  and  ./  **ro  o one 5 ante:  of  integration  to  he  evaluated  from 
tho  h  oondarv  c  oj  id  i.  t  i  our;  • 
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fa)  Structure  of  shook  wavo  without  react io 


With  the  assumption  that  the  substanco  Is  a  perfect  gas, 


P  -  pHT  =  EMT/u 


and  that  no  ohemloal  reaction  occurs, 


a  =  BL 


Equations  (46)  and  (47)  become 


Itu  4-RHT/u  -  3  s  y.  du/dy 


3^2  -fuj/li  -  -  7  -  jj  dl/dy 


Equations  (50)  ancl  (51)  aro  tho  basic  equations  for  a  pure  shock 
wavo  in  a  perfect  gas  —  such  a  wavo  as  might  be  initiated  and 
supported  by  a  moving  piston. 

Since  the  derivatives  on  the  right  sides  are  equal  to  aero 
before  the  detonation  wave  and  bohind  the  detonation  wave,  they 
may  be  set  equal  to  zero  and  the  two  equations  solved  simultaneous 
ly  for  the  temperature  and  velocity  before  and  after  the  detonatini 
wave.  This  Ib  tho  usual  procedure  of  tho  hydrodynamic -thermodjrrja^ir 
treatment.  On  tho  other  hand.  Equations  (50)  and  (51)  must  be 
solved  exactly  to  find  the  structure  of  the  wave  itself. 

?or  convenience  in  calculation,  we  mako  the  following 
changes  of  notation; 

Reduced  tomporaturo  3  =  f  HU2/J2 

Reduced  velocity  ^  =  u  li/3 

•Reduoed  distance  z  =  y  li/jj* 

£  =  Ii/2CV 

l»a  =  Z 

k  =  K>/f&y 


y 

r2,W,2 


L'he  following  differential  equations  are  obtained; 


d**>/dz  =  u)  +9/l*j  -  1 

ae/dz  =  l/k  £9  -  S(i-«)s  -  SaJ 
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If  Iflq.  (53)  1b  divided  by  i2q.  (52)  we  obtain  the  tempera ture- 
volocity  relation  for  a  particle 


d©  _  lAl  9  -  zLlA 

dto"  o)  +  0/iO  -1 


(54) 


:’or  the  spocial  oaso  k  =  1  +  2?,  aquations  (52)- (54)  oan  be 
solved  In  olosod  form,  for  If  0  Is  express oa  as  a  power  series  in 
o)  ,  the  series  fortuitously  terminates  (of.  Beoker).  Tho  result 

•£ 


0 


£ 

&S  +1 


(1  +  a  -<a2) 


(55) 


2^  1  to.ln  (H-fr>)  -  tin (*»-*>»«) 

z  =  J4-1  ,o  _i0. 


A  computation  of  the  temperature  and  density  within  a  shook  wave 
has  boon  carried  through  with  tho  aid  of  3qs.  (55)  and  (56)  and 
tho  assumed  numerical  values 

—  2. 

o  a  0.2  corresponding  to  -  6  cal  deg  mole  • 
k  =  1.4  which  is  about  correct  for  a  dilute  gas. 

a  2  0  which  means  that  the  detonation  velocity  is  very 
high  or  tho  initial  temperature  very  low,  since 

p 

for  low  rf0  a  =  2^1’q/D  approximately. 

The  results  of  the  computation  are  plotted  in  yig.  7. 


In  general  tho  conductivity  in  a  shock  vave  may  bo  expect e:l 
to  be  somewhat  greater  than  that  corresponding  to  k  5  1.4  •  It  is 
therefor©  desirable  to  solve  aqs.  (52)- (54)  for  any  value  of  k. 

An  attempt  to  solve  tho  equations  by  sotting  ©  equal  to  a  power 
series  in  <0  was  abandoned  becauso  the  series  did  not  converge  well. 
Beoourso  was  then  had  to  numerical  integration. 

The  derivative  in  Eq.  (54)  is  indeterminate  at  tho  boundary, 
*>inco  both  its  numerator  and  denominator  vanish.  However,  It  may 
bo  evaluated  by  tho  raothod  of  I* 'Hospital,  differentiating  both  with 
ve^poct  to  and  solving  tho  resulting  algebraic  equation  for 
c'G/d&e.  With  this  as  a  starting  point,  3q.  (64)  may  be  integrated 
.  umorically.  Tho  valuos  of  0  and  m>  so  obtained  permit  the  numerical 
integration  of  ISqs.  (52)  and  (53). 


OONFIDENW 

/■' 


L?i‘dTAKKE  CDDRT7fMATn  ^  Z 


S  jnio  oiirG  of  a 
a  pQi  foo'i;  rSilv  • 


shook  wav  a  without  ohemlcr.1  reaction  ..in 
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In  order  to  Investigate  the  of foot  of  inoroasod  thorma 
oonduotivlty,  tho  numerical  Integration  was  carried  through  for 
tho  numorloal  values 

S  Z  0.2 


6 

0 


and  tho  results  appoor  In  Jig.  6.  It  will  bo  seen  that  tho  only 
of feet  of  the  Increased  thormal  oonduotivlty  Is  to  spread  the 
toraporature-vvavo  farther  ahead  of  the  pressure-wave. 

In  order  to  Investigate  the  effect  of  a  higher  Initial 
temperature ,  numerical  integration  was  also  carried  through  for 
tho  values 

S  =  0.2 

k  =  6 
a  =  0.1 

and  tho  results  appoar  In  Jig.  9.  The  only  effect  of  the  higher 
initial  temporaturo  is  to  raiGo  tho  temperature -wave  slightly. 

Figures  7,  8,  and  9  are  all  drawn  to  the  same  scale  so 
they  are  directly  comparable.  In  oaoh  graph  temperature,  pressure, 
and  density  are  plotted  as  percentages  of  their  asymptotic  values. 
In  each  graph  the  unit  of  distance  £  is  about  45  Angstrom  units. 
These  shook  waves  are  comparatively"" wide  because  the  oaloulat lone  . 
wore  mado  for  a  gas;  for  a  solid,  the  waves  would  bo  much  narrower 

Discussion 


It  is  to  be  noted  that  in  each  of  tho  calculated  shook 
waves  tho  tempera ture-wavo  precedes  tho  pressure-wave.  This  appo; i? ?» 
to  be  a  general  phenomenon,  and  is  understandable  in  the  light  of 
tho  following  oonsldcrat ions • 

The  Impulses  which  form  tho  tomporaturo-wave  aro  transmit 
from  molooulo  to  rnolooulo  by  collisions.  Because  of  the  short 
distances  between  moloculos  in  tho  highly  compressed  state,  the 
fraotlon  of  the  distando  through  whioh  the  signal  must  be  oarried 
by  moving  molooulos  Is  small.  The  rate  of  transmission  of  the 
thermal  impulse  can  therefore  exceed  the  kinotio-theory  velocity 
of  tho  molecules  by  a  considerable  amount. 

The  donsity  (and  pressure  )-waves,  however,  are  built  up  by 
an  increased  oonoentrat ion  of  molooulos  themselves.  Tho  rate  of 
this  process  will  not  exceed  that  permitted  by  tho  kinet io-thoory 
velooity  of  the  molecules. 

If  tho  rato  of  transport  of  heat  always  exceeds  tho  rate  of 
transport  of  molecules,  we  may  safely  neglect  convective  transport 
jf  energy  in  the  detonation  wavo. 


CONFIDENTIAL 


FIGURE  & 


.-ssss  s.;  «• 

ina  unit  of  dlst.'.noo  z  la  0,IV^  ew** 


7. 


GONRDENTIM- 


Pl^T*MC£.  COORDINATE,  z 


structure  of  a  shook  wave  without  chemical  reaction  in 
a  perfect  gas.  Initial  temperature  higher  than  in  £i{?.  8. 
One  unit  of  distance  z  is  about  -V.  S. 
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(b )  ohook  wavo  with  chemical  reaction,  in  a  perfect  g, 
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Whon  a  ohomical  reaction  ie  occurring  within  tho  chock  wave 
the  basic  equations  for  tho  steady-stato  wavo  become 


i£u*f  RI.tT/u  -  J  =  yu  du/dy 


(57) 


0_2  -  (l-n)Q+uJ  4l  -  ■Jtt2  -  J  -  £  dj/dy 


i  (56) 


fjqs  •  (57)  and  (50)  differ  from  tho  oquations  for  tho  shock  wavo 
without  reaction  only  by  tho  inclusion  of  tho  term  (l-n)Q,  where 

&  is  the  heat  evolved  by  tho  reaction,  por  unit  m<.*ss, 
n  is  tho  fraction  of  matorial  left  unro acted. 


In  addition  thoro  will  bo  an  equation  for  tho  rate  of 
roaction,  which  v/ill  bo  for  a  first-order  reaction 


dn  . 
dt 


„  -A/i1 

2  o  ' 


(59) 


The  timo  derivative  in  Jflq.  (59)  can  be  ohangoti  to  the  space 
aorivatlvo  with  tho  aid  of  tho  relation  whioh  holds  for  the  steady 
stato 


dn 

dt 


■  u 


dn 

dy 


(60) 


For  convenience  in  oaloulation,  we  make  tho  same  changes  lr 
notation  as  before  (page  23),  with  tho  additional  onoe 

Q*  *  Q 
Z'  =  Zm/J 
A'  =  A  HH S/J2 


The  following  tliroo  difforontial  equations  oro  obtainod: 


d«i/dz  =  w  f  9/to  -  1  . ( 61 ) 

<16 /da  =  lA^S  “  Q'(l-n)  -  S  (X-tu>3  -SaJ  . (62) 

dn/da  =  -  -  Z'  o“A'/0  . (63) 
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liquations  (61),  (62),  (63)  must  bo  intogratod  simultaneous 
fhis  oan  bo  dono  by  numerical  integration. 
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The  result  of  the  integration  is  givon  in  Pig.  10  for  the 
conditions  ^ 

d  =  0.2  a  =  0.1  k  =  6 


Q*  =  0.1435  2*  =  10* 


A*  =  1.735 


which  correspond 


to  tho  boundary  conditions 


Z  =  2.5  1015  sec"1 
A  5  60,000/1.986  dog 
D  =  2.15  106  cm  sec"1 
Q  s  12.4  koal  molo-1 
Po  s  0.0012  gm  cm*3 
PQ  =  1  atmosphere 
=  4/3  s  2.27  10"4  poiso 


Probability  footor 
Heat  of  activation 
Dot on- t ion  velocity 
Heat  of  reaotlon 
Initial  density 
Initial  prossure 
Viscosity 


It  will  bo  seen  from  the  results  of  Pig.  10  that  a 
unimolocular  gas  reaction,  evon  with  a  high  activation  energy, 
will  go  to  completion  in  a  distance  less  than  that  required  for 
tho  systom  to  reach  its  highest  temperaturo.  Under  those  oonditl.c 
tho  rate-determining  stop  will  shift  to  some  other  process  such 
as  evaporation  or  activation  by  bimolooular  collision  for  gas 
reactions  (Cf.  disouscion  in  Port  I). 


(o)  Shock  wave  in  a  solid  oxplooivo,  with  evaporation  tho 
.rate-determining  step. 


When  a  shook  wavo  and  its  accompanying  chemical  reaction 
are  taking  place  in  a  solid  explosive,  the  equations  for  tho 
s toady-8 tat o  wavo  become 


u/V  -  M  ......( 64 ) 

Hu  ♦  P  -  J  s  f  \  du/d y  . . . . . .  ( 65 ) 

2  K> 

-  ( 1-n  )Q  -f-  uj /ll  -  *Ju  -  P  -  jt  d®/dy  ......(66) 
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The  pressure  P  1b  given  by  the  sum  of  tlio  prd&(4| 
the  gas  Initially  surrounding  the  grains  and  that  duo  to 
ovaporated  explosive: 


T/al 


r:pi^  +  a-n) 


V  -  (X 


whoro  is  tho  Initial  pressure  of  the  residual  gas, 

Is  the  Initial  spec  If  lv  volume, 

Ti  Is  the  Initial  temperature, 

n  Is  tho  fraction  of  explosive  remaining  in  the  solid  state , 
SL  is  the  ooVolume  (aesumod  constant). 

An  additional  equation  Is  roquirod  for  tho  rato  of  react iw., 
The  solid  is  supposed  to  be  grains  of  uniform  size,  and  It  is 
assumed  that  tho  reaction  rato  is  given  by  the  rate  of  evaporation 
of  moleculos  (or  portions  of  molecules)  from  tho  surface.  The  rato 
is  then  given  by 


-  dm/dt  -  g 


S  M 


whoro  m  is  tho  number  of  moluculos  per  gram  in  the  solid  state, 
g  is  tho  number  of  grains  per  gram, 

JT  is  tho  surface  area  of  a  grain, 

F  is  tho  area  occupied  by  a  molecule. 


Vpon  making  tho  substitutions 


m  -  nm 


there 


«  s  3T0/4‘fTro 
S  -  4<rrr^  n2/3 

b  =  (V>a)2/3 


m0  is  the  initial  number  of  moleculos  por  gram  in  the 

solid, 

VQ  is  the  initial  specific  volume  of  a  grain, 
rQ  is  the  initial  radius  of  a  grain, 
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the  rate  equation  booomoB; 

_  2/3  3  Vq/3  m  ^/R  -iHy'RT 

"  *»/**  8  n  57^7?  T  0  0 


'*%av7z 


(70) 


6ti 


whioh  may  bo  writ ton 

„  2/3  „ 

-  dn/dt  -  n  2  o 


(71) 


l>o  time  derivative  is  ohangod  to  a  spaoo  derivative  with  the  aid 
of  fcho  etoady-etate  rolation 


—  =  u  =  uv  M 

dt  dy  dy 


(72) 


Upon  making  tho  substitution  Vvs  V-Oc  ,  and  noglooting  £iSi£/5S.i. 
we  obtain  tho  three  different lj'.l  equations 


MyM.dV»/dy 

S  t,u2+  (l-n)Pi’A’  -  »* 

K  **/dy 

8  B^T  -  (l-n)Q  +  3*T»  -  ^V*2  -  3* 

-  dn/dy 

1  -AH%T  2/3 

-  t^t  zo 

o 


Jor  oonvonlenoe  in  calculation  wo  mako  changes  in  notation 
analogous  to  thoso  usod  before; 

Reduced  temporaturo  0  5  V  Wp/3  ,z 
Hoduood  velocity  W  s  V  M®/J  1 

Roduood  die  tone  o  z  -  y  U/j* 

S  =  R/2BV 

1  +  a  =  2H2i>,/J  ,2 

t  = 

Q1  =  Q  IH2/J,2(Tv 

z 1  =  a  f*/s  • 

A*  s  AH*  l^/j  ,2 

<X>  socjZ/J' 


onfidentim. 
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Upon  making  thoso  substitutions,  tho  three  differential  equations 
to  bo  Integrated  simultaneously  become 


d*/dz 
dO/dz 
-  dn/dz 


<a>  Ml-n)0//o-  1 


-  Q*(l-n) 


e 


-A*/© 


n3/3 


$  u-*>)2 


(76) 

(77) 
(76) 


While  analytic  solutions  have  not  been  obtained  for  the 
system  of  equations  (76),  (77),  (76)  in  the  general  oaso,  it  Is 
pocnible  to  solve  the  equations  oxaotly  in  the  front  of  the  wave 
rhere  the  amount  of  reaction  is  negligible.  The  solution  is 


<0 

Q 


z 


o 


+  02  o 


zA 


(79) 

(80) 


Solutions  over  tho  entire  rango  can  bo  obtained  for  specific 
substances  by  numerical  integration.  This  has  beon  carried  through 
for  TUT,  tho  results  of  tho  calculations  being  presented  in 
’Fig#  11.  Tho  numerical  integrations  are  somewhat  tedious  because 
tho  result  of  the  integration  is  sensitive  to  the  precise  value 
at  which  integration  is  star tou,  so  that  the  integration  must  be 
repeated  sovoral  timos  boforo  it  satisfies  all  tho  boundary  condi¬ 
tions.  Tho  final  rosult  in  Pig.  11  is  believed  to  bo  fairly 
acourato. 

Tho  numerical  values  assumed  for  TNT  in  this  calculation 

aro- 

Vx  S  0.629  om3  g"1 
D  -  6.9  103  om  8eo“* 

S  300°JC 

<X  =  0.4366  om3  g"1 
s  0.326  oal  g-1 

Mean  molecular  weight  of  products  =  25.9 
Q  s  770  oal  g"1 
Yg  =  0.6094  om3  g"1 
Tg  =  3293°K 

1^  s  2.39  10~*  oal  sec",L  om  (estimatod) 
ja.  ■  2  1(3®  poise  (estimated) 

QNF1DENTIAL 


1  I  4 


C  ioo7j 


-  36  - 


2ft 


FiGUI^E 


-4-0  -3o  -z© 

Pt^TAMf  e  Coordinate,  z. 


o 


Eltruoturo  of  a  shook:  wave  In  a  solid,  with  ohemioal  react  I?" 
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V„  =  0.6046  om3  g"1 

0  |  j  . 

kT/h  5  0.7553  10  *  boo"1  (estimated) 

ro  =  10  om  (ootloatod) 

AE*  =  11  koal  molo"1  (estimated) 


.  All  of  tho  proport ios  of  TNT  in  tho  abovo  list  except  those 

marked  nest imatod”  wore  taken  from  a  calculation  of  the  detonation 
properties  of  TUT  by  Cook.  The  only  one  of  tho  estimated 
proportion  v/hioh  roquiros  comment  is  the  grain-radius  r^  which 

was  taken  as  10~7  cm;  this  value  is  surely  an  extreme  lNtoor  limit , 
aad  since  tho  rate  of  reaction  Is  inversely  proportional  to  the 
initial  grain  else,  this  calculation  roprosents  an  extreme  upper 
jimit  for  tho  rate  of  reaction  by  tho  evaporation  moohanism. 

Ono  unit  of  z  in  ?lg.  11  is  0.2  ingetrjrin  units,  which 
would  make  tho  roaoTlon  zono  (and  temperature-wave)  about  20  or  30 

Sng.itrffm  units  long  in  all. 

Discussion 


Several  features  of  these  rosults  doserve  comment; 

(1)  Tho  wave-front  of  tho  temperature  wave  is  comparatively  broad 
and  of  a  gradual  slope,  tho  temperature  appearing  to  riBO  as  the 
reaction  proceeds.  See  the  discussion  on  page  21. 

(2)  The  pressuro-wavo,  on  tho  contrary,  is  narrow  and  steep.  Just 
as  in  other  shook  waves,  tho  pressure  front  lags  far  behind  tho 
tomperaturo  front.  The  effect  of  pressure  changes  on  the  reaction 
is  negligible,  and  tho  pressure  front  is  too  far  behind  tho  temp- 
orature  front  to  have  any  important  effoct  in  initiation  of  tho 
reaction.  When  the  pressure  front  arrives,  it  does  produce  a 
sudden  jump  in  tho  tomperaturo. 

(3)  It  has  boon  proposed  that  detonation  roaotions  should  have  a 

wave  form  like  that  in  ?ig.  12,  the 
temperature  and  pressure  passing 
through  an  exceedingly  high  maximum 
just  behind  the  front  of  the  wave,  and 
approuohing  their  asymptotio  values 
from  above.  The  differential  equation? 
do  Indicate  tho  existence  of  such  a 
maximum.  However,  the  calculations 
which  Jig.  11  illustrates  show*  that 
tho  maximum  in  pressure  Is  not  large, 
while  the  maximum  in  temperature  is 
imperceptible  on  tho  graph.  These 
maxima  are  therefore  probably  not  of 
sufficient  magnitude  to  be  of  c.ny  groat 
importance . 
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If  explosive  charges  ore  of  suffLoient  diameter  that  the 
wave  does  resemble  olosely  a  one-dimensional  wave,  It  Is  possible 
to  obtain  detonation  velocities  very  olose  to  those  theoretically 
calculated*  However,  actual  explosive  charges  are  often  of  small 
enough  diameter  that  a  considerable  port  of  the  energy  of  the 
detonation  wave  is  dissipated  laterally  —  producing  slowing  or 
eventually  stopping  of  the  detonation  wave.  The  diameter  at  whioh 
this  offeot  becomes  appreciable  may  be  of  the  order  of  inohes  for 
now-energy  explosives. 

It  may  be  supposed  a  priori  that  the  detonation  wave  ceases 
because  the  chemical  react  Ton  which  was  its  support  has  ceased. 

At  a  point  where  there  is  no  detonation,  there  is  no  longer  any 
?hemioal  reaction. 


As  the  discussion  of  Part  I  indicates,  the  pressure  has 
little  or  no  effect  on  the  rate  of  reaction.  We  therefore  suppose 
the  ohemloal  reaction  becomes  excessively  slow  only  because 
'*•1  ?  temperature  has  fallen  excessively  low.  When  the  explosive  is 

mated  with  tho  aid  of  a  booster,  the  temperature  at  the 
•.equator  end  is  sufficiently  high  (if  the  booster  is  supposed  to 
be  a  high-energy  explosive)  that  the  reaction  will  go  in  the  solid 
phase.  But  as  the  detonation  wave  proceeds,  if  heat  is  lost  rad¬ 
ially  the  back  part  of  the  detonation  wave  will  be  oooled.  Phis 
will  oause  the  reaction  to  beoome  slower  in  this  region,  whioh 
will  then  cool  the  region  immediately  preceding  it.  The  result 
will  be  that  a  region  of  low  pressure  and  temperature  (rarefaction 
wave)  will  spread  forward  until  it  reaches  the  front  of  the 
detonation  wave  and  the  detonation  will  fail. 

The  effect  of  temperature-lowering  on  the  reaction  rate  will 
be  particularly  noticeable  if  the  reaction  has  to  go  homogeneously 
if  at  all,  because  it  has  no  surface  to  permit  the  heterogeneous 
reaction.  This  will  be  true  for  cast  explosives,  and  to  a  lesser 
extent  true  for  explosives  composed  of  large-sized  particles. 

The  high  activation  energy  of  the  homogeneous  reaction  will  oause 
the  reaction  to  come  to  a  complete  standstill  even  at  moderately 
high  temperatures  (of.  Pig.  4). 

Tho  complete  problem  of  the  rate  of  dissipative  loss  in  a 

wave  with  chemical  reaction  is  a  difficult  one  which  has  not  been 
solved.  However,  the  magnitude  of  dissipative  lose  as  a  function 
of  the  amount  of  radial  expansion  oan  be  studied.  There  are  at 
least  three  possible  ways  in  which  expansion  might  produce  a 
lowering  of  temperature: 


If  the  chemical  reaction  in  a  particular  region  has  only 
gone  partially  to  completion  when  the  explosive  has  expanded,  the 
remaining  reaction  will  take  place  as  if  the  initial  density  of 
loading  of  the  explosive  were  less.  This  might  effect  both  the 
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pressure  and  the  temperature  whioh  the  oomplete  reaction  would 
produoe. 

The  experimental  fact  is  that  while  the  initial  density  of 
loading  has  a  marked  effect  on  the  detonation  velocity  (Fig.  13 
gives  the  effect  on  detonation  velocity  for  some  common  explosives) 
and  on  the  detonation  pressure  (Pig.  14.1/  gives  the  effect  on 
pressure ),  it  has  substantially  no  effeot  on  the  detonation 
temperature  (?i g.  16). 

This  result  should  have  been  expeoted,  for  the  detonation 
temperature  is  determined  primarily  by  the  heat  evolved  by  the 
reaction  and  the  heat  capacity  of  the  products,  neither  of  whioh 
is  particularly  sensitive  to  density.  Uore  exactly,  the  detonation 
temperature  for  a  solid  explosive  is  given  by 

+■  Q/&,r 


-  d*/AV2 ) 


w'iwre  H  is  the  number  of  molecules  produoed  from  a  single  molecule 
i" .7  hc~oxplosion,  and  tho  other  symbols  have  their  usual  meaning 
It  will  bo  noted  that  the  density  appears  explicitly  only 
I  -he  derivative  of  the  oovolume.  The  effeot  of  density  would 
'jrxor&fore  be  small. 

(2)  Effect  of  adiabatic  expansion 

If  tho  reacted  (or  partially  reacted)  explosive  behind  the 
detonation  front  expands  so  qulokly  that  heat  loss  to  the 
surroundings  may  be  neglected,  we  can  calculate  the  temperature 
drop  caused  by  the  adiabatic  expansion. 

For  a  substance  which  obeys  the  imperfect  gas  law 


P  (V  -0O  =  2T 


with  the  covolume  constant,  the  adiabatlo  law  is 


2?(V-fc)CP/Cv  =  oonbt, 

r/c_ 

T  (V  -  OO  v  =  oonst , 


from  whioh  the  temperature  and  pressure  after  ^ny  given  expansion 
may  be  computed.  ThiB  has  been  done,  with  the  assumption  that 
,  the  oovolume  o<  is  0.85  of  the  detonation  volume  Vo.  The  results 
are  plotted  in  Fig.  16. 
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(5)  Effoot  of  Inertia  of  the  pace* 


If  work  1b  done  against  the  inert  La  forces  due  to  the 
oaso  surrounding  the  explosive  charge,  energy  will  be  dissipated 
thereby  and  the  temperature  will  fall. 

The  results  for  this  process  have  been  oaloulated  by 
numerical  integration  by  5.1.  Taylor,  for  a  oylindrioal  explosive 
burrounded  by  a  heavy  oase,  and  the  case  does  not  begin  to  expand 
until  the  detonation  wave  has  passed  and  the  reaotion  is  over. 

The  temperature  and  pressuro  oaloulated  in  this  way  for 
-arious  degrees  of  expansion  are  plotted  in  Pig.  17,  the  data  for 
which  were  taken  from  Table  I  of  Taylor’s  report.  The  temperature 
■i.as  not  given  in  Taylor’s  Table  I,  but  was  oomputod  from  his 
■values  of  prossure  and  density  with  the  aid  of  the  equation  of 
*?tate. 

The  calculations  of  Taylor  on  which  Pig.  17  is  based  include 
>*utla  the  effect  of  inertia  and  of  adiabatio  expansion. 


Discussion 


It  is  apparent  that  the  offsets  of  radial  expansion  will 
lower  the  temporaturo  bohind  the  front  of  a  detonation  wave  — 
tAnd  the  magnitude  of  tho  temporaturo  effect  is  l«..rgo  enough  to 
have  a  major  effect  on  tho  rate  of  reaction. 

Por  oxample,  an  expansion  to  twice  tho  initial  volume  will 
lower  the  temporaturo  to  about  of  its  initial  value  (Pig.  17). 
If  the  initial  temperature  was  4000°K,  this  cooling  would  be 
sufficient  to  bring  it  to  800°E.f  which  would  effectively  stop 
not  only  tho  homogonoous  (solid)  roaction  tut  also  the  hetero¬ 
geneous  (ovaporation)  roaction  and  so  bring  all  decomposition 
of  tho  explosive  to  a  standstill. 
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PART  IV.  THIS  DETONATION  WAV  IS  IN  TILS  NON-STEADY  STAT 


Having  doterminod  the  distribution  of  prossuro,  temperature, 
density  and  volooity  in  tho  steady-state  dotonation  wave,  wo  may 
*  vlsh  to  dotermlno  tho  same  properties  for  the  non-steady  state. 

Jor  example,  wo  may  wish  to  hr.vo  tho  wave  pass  from  ono  medium  into 
another  with  different  properties  and  follow  its  progress  there. 

«  In  principle,  it  is  possible  to  do  this  by  numerical 

.r; ^ration  of  tho  general  equations  of  conservation,  introducing 
.ihosen  conditions  dlscontlnuously  at  some  point  on  tho  spaco 
ocoriinato.  Any  given  distribution  may  be  differentiated  graphically 
or  tabularly  to  obtain  all  the  nooessary  terms  in  the  equations 

■&*"s(eu)  . (»i 

. ,«,) 

H  ?  **  '  r  W 

21.  —  fp-u^AM  K>  "br  _£Jn  . (06) 
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rich  tho  usual  auxiliary  equations 


P  *  f(p  %  T)  (Equation  of  state)  . (8*7) 

E  =  CyT  -  (l-n)Q  . (88) 

dn/dt  =  ;  P(nt  T)  (Bate  law)  . .(89) 


Having  obtained  tho  time  derivatives,  it  would  bo  possible  to  find 
now  values  for  tho  variables  by  relations  such  as 

2>  =  T+  eto.  . (90) 


Repetition  of  this  process  will  then  show  graphically  the  effect 
on  a  detonation  wave  of  any  sot  of  conditions  we  choose  to  impose. 
Tho  wave  oould  bo  followod  from  ono  explosive  into  mother,  as 
for  a  booster.  The  method  oould  also  follow  the  thermal  initiation 
of  an  explosivo  by  a  givon  temperature  distribution. 

Tho  mothod  has  boon  tried  by  the  writers  in  a  preliminary 
way,  but  tho  caloulation  involved  is  oxcossive  and  no  final  results 
have  boon  obtainod  in  this  way.  An  attempt  to  evaluate  the  deriva  ¬ 
tives  by  oxpressing  T  as  a  power  series  in  t_  (see  Informal  Progress 
Report  Z)  was  abandoned  bocauso  it  did  not  converge  fast  enough  to 
permit  any  usoful  calculations . 
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